Dynamics of colloidal particles in ice 
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Abstract 

We use X-ray Photon Correlation Spectroscopy (XPCS) to probe the dynamics of colloidal 
particles in polycrystalline ice. During freezing, the dendritic ice morphology and rejection of 
particles from the ice created regions of high-particle-density, where some of the colloids were forced 
into contact and formed disordered aggregates. We find that the particles in these high density 
regions underwent ballistic motion coupled with both stretched and compressed exponential decays 
of the intensity autocorrelation function, and that the particles' characteristic velocity increased 
with temperature. We explain this behavior in terms of ice grain boundary migration. 

PACS numbers: 82.70.Dd,64.75.Xc 
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Solidification of the solvent phase of a colloidal suspension occurs in a wide variety of 
natural and technological settings. Owing to the rapidly expanding domain of materials 
applications, and relatively simple and inexpensive processing methods, a variety of solidifi- 
cation processes are receiving intense theoretical and experimental study Under a wide 
range of conditions, as the solvent freezes the interface between the liquid and solid solvent 
ramifies, guiding the particles into a variety of macroscopic morphologies. At the lowest 
freezing velocities all of the particles are rejected and pushed ahead of a macroscopically 
planar solid-liquid interface, while at higher velocities a dendritic solid-liquid interface aligns 
the particles into microporous structures, or solid lenses segregate the particles into peri- 
odic layers 2, j|. At the highest growth rates the solid engulfs the particles individually. 
These various regimes have many applications including purifying water [4|, creating tissue 
scaffolds {5], understanding ground freezing jg], and building composite materials The 
progression through the various regimes and their properties are fundamentally important, 
but involve a hierarchy of poorly understood cooperative phenomena. One of the factors 
limiting theoretical understanding is the dearth of experimental studies exploring the physi- 
cal processes at work on the particle-scale. Here, we provide the first investigation of particle 
dynamics within frozen colloidal suspensions. 

We used X-ray Photon Correlation Spectroscopy (XPCS) to examine the dynamic be- 
havior of spherical silica colloids in ice. The scattering reveals ballistic dynamics combined 
with a non-exponential decay of the intensity autocorrelation function (ACF). This combi- 
nation with a compressed exponential decay is commonly observed in light scattering from 
soft materials. Uniquely, in our experiments the decay is slower than exponential at small 
scattering vectors and faster than exponential at large scattering vectors. While ballistic 
dynamics combined with a stretched exponential have been reported previously in one other 
system, ours is evidently the first observation of a transition from stretched to compressed 
exponential behavior with increasing scattering vector. As such, it offers insight into the 
source of ballistic particle dynamics and non-exponential decay of the intensity ACF in 
non-equilibrium colloidal materials. 

Samples of colloidal silica spheres (Polysciences, Inc.) dispersed in deionized water were 
solidified, and the composite material interrogated via X-ray scattering. The particle radius 
R = 32 nm, polydispersity z = 18%, and initial unfrozen particle volume fraction 2% 
were determined from small angle X-ray scattering (SAXS) Q]. The samples were contained 



in an approximately 400 /im thick, temperature-controlled sample chamber that produced 
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Observations were made at several 



a cylindrical isothermal region 2 mm in diameter 
locations within the isothermal region of the samples. 

The solvent was frozen by lowering the temperature at about l°C/s. Freezing usually 
occurred near T = — 25°C Therefore, the water in all samples was highly supercooled when 
ice growth began, resulting in an unstable solidification front and a cellular or dendritic 
ice growth morphology 3, 9, 10]. Indeed, when we observed samples freezing under similar 



conditions with video microscopy the ice growth was dendritic [8j. As expected 10], the 
colloidal particles were rejected into the regions between the relatively pure ice dendrites 
resulting in mm-scale linear regions of high particle density tens of /im wide separated by 
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regions of low particle density of approximately the same width 

After freezing, we increased T in steps as small as 0.05°C from T = — 2°C to T = 0°C, and 
obtained sets of images of the scattered X-ray intensity / (q, t) at various scattering vectors 
q and times t. We performed the X-ray scattering at sector 8-ID-I of the Advanced Photon 
Source ll]. The beam size was 20 /im by 20 fim, ensuring that the illumnated volume of 
the sample contained many millions of particles. We used the multispeckle XPCS method 
to obtain the intensity autocorrelation function g 2 (q, r) = (I (q, t) I (q,t + r))/ (I {q)) 2 from 
the azimuthally-averaged intensity I(q,t). Here (• • •) represents an ensemble average over 
the detector pixels and r is the delay time between the two frames for which the correlation 
is calculated. The intensity ACF contains information about the particle dynamics through 
the shape and rate of its decay [12] . The length-scale of the dynamics probed is related to q; 
by using X-rays, we access length-scales comparable to the size of the individual particles. 

Figure [1] shows several ACFs at q ~ 0.1 nm _1 (near the peak in the static scattered 
intensity) for a sample at T = -0.60°C, T = -0.50°C, T = -0.40°C, and T = -0.30°C. 
These data represent the typical forms of g 2 (q, r) observed: either a single decay, or two 
decays in which the second is extremely stretched and of small amplitude. For data with 
two decays, we are only concerned with the decay at small r. As the temperature increases, 
the decay rate also increases, so that the decay time becomes shorter. 

In order to quantify the changes in decay time and shape, we fit each g 2 (q, r) to 



a Kohlrausch- Williams- Watts (KWW) expression 13] or a combination of two KWW 
expressions. Single decays were fit with g 2 (q, r) = b + Cexp [— 2 (Tr) a ], where b is 
the baseline, C is the contrast, T is the decay rate, and a is the stretching (< 1) 
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FIG. 1. Adjusted intensity ACFs [<?2(<Z> T ) — b]/C versus delay time r from sample 2 (of four) 
at T = -0.60°C (circles), T = -0.50°C (squares), T = -0.40°C (triangles), and T = -0.30°C 
(diamonds) with the corresponding stretched or compressed exponential fits (solid curves). 



or compressing (> 1) exponent. Similarly, we fit double decays with (q, r) = 1 + 
C {(1 — (3) exp [— (T F r) aF } + /3exp [— (T s r) as }} 2 , where C, T and a are as above. The 
"partition coefficient" (3 describes the relative strength of the two exponential decays, and 
the subscripts denote the first (F) or second (S) decay. As illustrated by the solid curves 
in Fig. [1], most data were fit well with one of these functions, but some were too noisy or 
had features that did not fit into this analysis framework. We ascribe this variability to 
the spatial heterogeneity produced by the inherently stochastic nature of the ice nucleation 
process, the unstable ice growth morphology, and the process of ice crystal coarsening in 
the polycrystalline ice. 

Both T and a resulting from these fits vary with q. Figured shows examples of T versus 
q for a sample at T = — 1.60°C, T = — 1.10°C, and T = — 0.70°C, all of which increase 
linearly with q. Figure [2b shows examples of a versus q from the same data. In all cases, 
a is between 0.5 and 1 at low q and increases to between 1 and about 1.5 at higher q, 
leveling off for q ^ 0.12 nm _1 . For comparison, the measured structure factor from SAXS 
has a peak at q = 0.123 nm" 1 and the particle size corresponds to q — tt/32 nm = 0.0981 
nm" 1 [8]. This behavior indicates ballistic particle dynamics and can be described with a 
distribution of particle velocities similar to a Levy stable distribution, which has a power 
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FIG. 2. In (a), examples of decay rate T are plotted versus scattering vector q from sample 2 at 
T = — 1.60°C (circles), T = — 1.10°C (squares), and T = — 0.70°C (triangles) with their respective 
linear fits (solid lines). In (b), the corresponding exponents a are plotted versus q. 



law tail [14]. The characteristic velocity v c of the particles is the slope of a linear fit to V 
versus q (solid lines in Fig. |2^l). We have fit T from the data analyzed in this framework to 
find that v c generally increases with increasing temperature as shown in Fig. |3j The orders 
of magnitude difference between the two curves can be ascribed to the variability introduced 
by the freezing process. 

While ballistic dynamics combined with a compressed exponential decay of #2 (<?, t) 



has become a common observation 14j4l7l]. the combination of ballistic dynamics with a 
stretched exponential decay is rare 3] and a transition from stretched to compressed ex- 
ponential decay with increasing scattering vector has not been reported previously. Either 
stretched or compressed decays could result from a continuous time random walk in which 
the size of the steps performed by the scatterers is distributed according to the Levy alpha- 
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FIG. 3. The characteristic velocity v c versus temperature T for sample 2 (circles) and sample 3 
(squares) . 



stable distribution. This explanation is not specific to a particular material or system, but 
ascribing a physical origin to the neccessary step size distribution can be difficult. For col- 
loidal gels, considering the relaxation of stress dipoles induced during disturbances to the 
material (e.g . shear or loading into a container) produces the values of a obtained experi- 
mentally 16, l^. However, few other systems exhibiting this type of behavior have such a 
complete description. 

Here we suggest that the high-particle-density regions, which dominate the scattering, 
contained colloidal aggregates that formed when the solutions froze and that grain boundary 
motion arising from coarsening of the polycrystalline ice exerts forces on these aggregates 
resulting in the observed dynamic scattering. From SAXS performed simultaneously with 
XPCS, we know that the particles in the high density regions were close-packed and touching 
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8j. The forces associated with particle rejection from the ice front during freezing are 
sufficient to overcome the electrostatic repulsion between the particles and bring them into 
contact [20]. Once in contact, the van der Waal's attraction between the particles caused 
them to bind to each other forming an amorphous colloidal solid. The attractive energy for 
these particles is estimated to be 9 x 10~ 20 J, or about at T = 273 K, so thermal 

energy alone would be insufficient to completely disperse the aggregates on experimental 
time-scales. Indeed, using video microscopy we observed high-particle-density aggregates 
up to hundreds /xm in size sedimenting through the melted samples. SAXS measurements 
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indicated that about half of the particles initially in the solution sedimented [§J. Thus, the 
grain boundaries between ice crystals likely contained a mixture of aggregates and single 
particles. 

In order for the polycrystalline ice to coarsen, the aggregates and particles in the grain 
boundaries must deform under stresses transmitted from the ice crystals. We estimate the 
strength of these forces from the pressure AP driving the ice crystal coarsening: AP = 'jgbH, 
where jgb is the surface free energy of the grain boundary and k is the curvature of the 



boundary [21]. Taking 7^ = 0.029 J/m 2 , which is the value for an ice-water interface 22], 



and k = 1/R C = 10 4 m _1 , where R c = 100 /im is an estimate of the ice grain boundary radius 
of curvature 2jJ, gives AP = 290 N/m 2 , which is distributed across the particles adjacent 
to the ice. Assuming that the pressure is distributed across about half the surface of a 10 /im 
radius aggregate, the total force on the aggregate is about 3.6 x 10~ 7 N. At = 0.6, this 
amounts to a force per particle on the surface F p of about 2 x 10~ 12 N. In order for the 
aggregate to deform, we assume that a particle must move a distance equal to one particle 
diameter (64 nm). The work performed by F p in moving a particle this distance is about 
1 x 10 -19 J, slightly greater than the van der Waal's energy binding the particles together 
(9 x 10~ 20 J). Therefore, grain boundary motion induced by coarsening of the polycrystalline 
ice can deform the colloidal aggregates. 

Although other stress sources are possible, grain boundary motion from coarsening pro- 
duces characteristic velocities of the correct order of magnitude. The grain boundary velocity 
is given by v g b = M^gbK, where M is the grain boundary mobility. The product M^ g b for 



ice at T = — 5°C ranges from about 10~ n to 10~ 14 m 2 /s 24]. Thus, for the value of k 
given above, v g b could be between 0.1 and 100 nm/s encompassing the observed range of v c . 
Variations in grain size and mobility among the samples can account for the factor of 100 
difference between the characteristic velocities shown in Fig. [3J Furthermore, grain bound- 
ary migration can account for the increase in the characteristic velocity with temperature 
through the Arrhenius behavior of the grain boundary mobility 24J]. 

Determining how the stresses produced by grain boundary migration would affect the 
motion of particles in the grain boundaries is challenging; we discuss two complementary 
scenarios. It is unlikely that the stress resembles the dipolar stress sources proposed to 
explain compressed exponential decay and ballistic motion [3]. Instead, the very high 
particle density may cause the colloids to behave more like a granular material (e.g., 
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so that forces are distributed along force chains [26] . Such a direct transmission of the grain 
boundary motion would produce ballistic particle motion at about the same rate as the 
grain boundary was moving. However, the distribution of the force would cause particles to 
experience stress only sporadically At larger length scales (low q), there would be a wider 
range of decay times (smaller a) because the sporadic stress may be greater or less than that 
required to cause the necessary particle displacement for decorrelation of the intensity ACF. 
Whereas, at smaller length scales (large q) almost any stress would deform the particles 
sufficiently for decorrelation and thus there would be a narrower distribution of decay times 
(larger a). However, the stress is still imposed sporadically leading to a non- negligible 
possibility of very long waiting times between imposition of sufficient stress, and hence a 
power law form. 

We can connect this type of behavior with prior theoretical work on solidifying colloidal 
suspensions and self-filtration. In colloidal suspensions near random close packing, the os- 
motic pressure, and hence the generalized form of the Stokes-Einstein (mutual) diffusivity, 
diverge as — > <fi rcp Then, the character of diffusive relaxation of particle con- 

centration gradients becomes extremely long ranged such that very near the divergent limit 
diffusion becomes effectively "instantaneous." Thus, imposing a force on the close-packed 
colloids in the grain boundary leads to immediate and long-ranged particle motion, i.e. force 
chains. Such behavior of the diffusivity can be understood to drive ballistic motion. 

Given the increasing study of high density and non-equilibrium systems, more examples 
of ballistic dynamics coupled with non-exponential decays are sure to emerge. Many of these 
will not fit into the developed framework and the manner of force transmission should be 
considered. For example, a similar type of stress distribution along force chains may explain 
the ballistic dynamics combined with a stretched exponential decay observed in particle- 
stabilized emulsions 18j. Our experiments highlight the need for more general studies to 
quantify how different types of stresses manifest themselves in particle motion and scattering. 
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